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What is measured in the scanning gate microscopy of a quantum point contact? 
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The conductance change due to a local perturbation in a phase-coherent nanostructure is cal- 
culated. The general expressions to first and second order in the perturbation are applied to the 
scanning gate microscopy of a two-dimensional electron gas containing a quantum point contact. 
The first-order correction depends on two scattering states with electrons incoming from opposite 
leads and is suppressed on a conductance plateau; it is significant in the step regions. On the 
plateaus, the dominant second-order term likewise depends on scattering states incoming from both 
sides. It is always negative, exhibits fringes, and has a spatial decay consistent with experiments. 
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Scanning gate microscopy (SGM) has been intensively 
used during the past decade to investigate electronic 
transport in nanostructured two-dimensional electron 
gases such as quantum point contacts [IHl] , quantum bil- 
liards iS], and rings [HIE]- In this technique, a charged 
tip that locally influences a device's electrons is scanned 
over the surface; i.e. the perturbed nanostructure's con- 
ductance is a function of the tip position. The SGM tech- 
nique is of great fundamental interest since it provides a 
wealth of data that depends on the microscopic electron 
transport. It can generate very detailed sample charac- 
terization and more precise information of the disorder 
configuration than a traditional transport measurement. 

Of importance in SGM studies is the precise physical 
interpretation of data. Are the measurements sensitive 
to the electron flow as often just stated or, perhaps, the 
local electron density? Whereas experiments and numeri- 
cal simulations have yielded conductance change patterns 
caused by local perturbations closely related to calcu- 
lated electron flows or local electron densities [H [31 |S1- 
110] , the precise relationships among these quantities and 
interpretations remain far from obvious. Moreover, it is 
recognized that in a linear response framework the cur- 
rent density is not uniquely defined, and its study merely 
provides a way to visualize the structure of the scatter- 
ing states at the Fermi energy [11. Thus, our goal is to 
provide a clear-cut relationship between the conductance 
change in a SGM setup and the unperturbed scattering 
states which is valid for an arbitrary coherent structure. 

Gonsider a system with a Hamiltonian H = Hq + V, 
where Hq represents the unperturbed structure and V 
the tip potential. Though quite general, the main focus 
is on a quantum point contact (QPG). This paradigm ex- 
hibiting conductance quantization [12H16| has been well 
studied by SGM [HIJ- In addition, recent work suggests 
the relevance of this powerful technique for probing the 
nonlocality of electronic interactions P^THTO] . 

The scattering theory of quantum conductance [2Ql EI] 
assumes that the leads are disorder-free, confined in the 



transverse (y) direction, and semi-infinite in the longi- 
tudinal (x) direction. The electron states in the leads 
are products of quantized transverse wavefunctions 4>a{y) 
labeled by the index a with transverse energy £a and 
plane waves propagating in the x direction with wave- 
vector magnitudes ka- The total electron energies are 



-M 



h k^/2M, with M the effective electron mass. 
The N transverse momenta which satisfy this relation- 
ship with fc^ > define the 2N propagating channels of 
the leads with energy e. Denoting r = {x,y), the incom- 
ing (-) and outgoing {+) lead states are defined as 
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For x > 0, the index 1-^2 and the opposite sign is 
taken in the argument of the exponential. An infinitesi- 
mal negative (positive) imaginary part is given to ka for 
incoming (outgoing) lead states. For simplicity, consider 
a confining potential that is x-independent in the asymp- 
totic region even though any separable potential can be 
treated if allowance is made for an x dependence of ka- 

The scattering states corresponding to an electron in- 
coming from the left lead 1 (right lead 2) with energy e 
in the mode a are given in the asymptotic regions by 
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Then the 2Nx2N scattering matrix S, relating incoming 
and outgoing fluxes, can be written in terms of the NxN 
reflection and transmission matrices from lead 1 (2) as 



r t' 
t r' 



(3) 



The chosen lead state normalization corresponds to an 
incoming flux e/h and ensures the orthonormality of the 



scattering states [201 HI] ■ The zero-temperature conduc- 
tance of the unperturbed structure in units of the con- 
ductance quantum 2e^//i is given by g^^' = Tr[i'''i], where 
Tr denotes the trace over the modes. 

The Lippmann-Schwinger equation for the perturbed 
wave function Xi ^ ai'"^) using the retarded Green function 
tj(°) associated to Hq is 



^l,e,a 



W = ^tlW + / dr g^"Hr,r,e)V{r)x[^Jjr) 



(4) 



It is consistent within Hnear response theory to calcu- 
late the conductance change beginning with the current 
change carried by "^\J^{r) with the spectral decompo- 
sition of t/('') in the scattering wave-function basis. To 
first order in V, the integration over a cross section Sx 
on the scatterer's right [22^ gives 
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where the matrix element of the perturbing potential in the scattering state basis is 

VkUe-,e) ^ jdv ^\l>{v) V{v) M/^XUr) • 
In addition, a shorthand notation is introduced for the following quantities involving the unperturbed states: 

N 
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(7) 



N 



t*aa exp [-i(A:a + k^ )x] + ^ 



&=i 




^ ' C''fcaexp[i(fc^-fcf, )x] 



with iV = min{iV,iV}. 

The £ integral of Eq. (Is]) has a principal part and a S- 
function contribution. Since the current correction has to 
be independent of x, the integral can be evaluated in the 
limit a; — >■ oo. The principal part contribution is simply 
equal to that of the (5-function contribution. The result 
is thus for £ = £, and forthwith the energy arguments are 
dropped. The current change associated with the mode 
a at energy e is given by 



Ii:U-T.^^{i*'t^'''--'i'^'^')a.a}- (8) 
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The first term vanishes when summed over a. Within 
linear response in the bias voltage, the first-order change 
in the zero-temperature conductance is 
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(9) 



This result is valid for the general situation where quan- 
tum transport through a scatterer is modified by a weak 
perturbation [53|. The matrix rH' depends only on the 
unperturbed scatterer, while the tip's effect appears in 
the V matrix elements. The conductance change is not 



simply given by the current density or charge density at 
the tip position but is proportional to products involving 
a matrix element of V with two scattering states, one cor- 
responding to an incoming electron from the left and the 
other from the right lead. The conductance's sensitivity 
to electrostatic potential variations has been considered 
in a one-dimensional geometry in Ref. [24] . The example 
of a (5-function barrier perturbed by a local tip can be 
analytically calculated and agrees with Eq. ^. 



For a QPC, rH' is appreciable only in the vicinity of 
the conductance steps and is suppressed on the plateaus. 
A basis transformation into the eigenmodes of tH singles 
out the linear combinations corresponding to the modes 
that propagate through the constriction [T^. Only at 
these channels' openings is rH' significant. The first- 
order conductance correction [Eq. ^] is thus the rele- 
vant term in the step regions, but to capture the plateau 
regions' dominant correction |25j, the second-order con- 
tribution to the scattering wave functions in Eq. Q must 
be taken into account. Continuing with the same method 



gives the second-order conductance correction 

g(2) ^ _4^2 r^^ [^ti Vl.2y24 ^ ^,^, -^2,1-^1,2 

+ Rc{rV(v2>2-^2,i_-^;24via)|] (^q) 

This result, together with Eq. Q, provides the answer 
to the question posed in the title. 

Like g^^\ g^^^ depends on the scattering amplitudes of 
the QPC and the scattering states incoming from both 



sides. In contrast to g 



(1)_ g{2) 



contains nonvanishing 



terms with pref actors tH and t'H' . Since rH' is sup- 
pressed on the plateaus, g'^^^ is expected to be the dom- 
inant correction there. Furthermore, the second-order 
terms contain contributions from all effective transmit- 
ting channels, in contrast to g^^\ which is significant only 
in the steps and dominated by the opening mode. Thus, 
the transverse shape of the channels should be observ- 
able directly with SGM on the steps without subtracting 
the lower mode results as done in the plateau data anal- 
ysis [T]. 

Crucial insights into conductance quantization of a 
QPC have been gained through the analysis of simple 
models, such as the saddle potential [TB] for the smooth 
electrostatic potential of the constriction. In addition, 
describing the saddle potential by a double harmonic os- 
cillator V{r) = Vo + ^{ujyy'^ — ^x^'^) allows for an exact 
solution |26j . The transmission and reflection amplitudes 
are diagonal in the mode indices with 
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where a{£) — £ -|-argr(i + i£) — £\n\£\ in terms of the 
dimensionless energy £a — 2[e — fiujy{a + 5) — Vo]/TiU!x- 
In the limit ujy 3> Wj;, a good conductance quantization 
is achieved [16]. As expected for an arbitrary QPC, \r*t\ 
is appreciable only at the steps, and on the A^*'^ plateau 
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Assuming a local tip with V{r) = VT<^(r — i"o) and using 
the semiclassical form of the scattering wavefunctions, we 
have 
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for the propagating modes, where I — y/h/2Muj^ and 
/C = ^y£a + (xo/^iy. The a*^ eigenstate (f>a (y) of a one- 
dimensional harmonic oscillator with frequency cuy yields 
the lobes in the transverse direction that are ubiquitous 
in experiments and simulations. V^'^ has real and imagi- 
nary parts that oscillate as a function of xq with an over- 
all decay proportional to Xq far from the constriction. 
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FIG. 1: Top: dg/dVr of the conductance change with tip 
voltage for a hard wall QPC and a (5-tip at the first conduc- 
tance step (point s of the conductance curve shown in the 
inset) as a function of the tip position xq. The solid line is 
for yo — (in the center); the dashed and dotted hues are 
for yo — 4: and 8, respectively. Inset: The conductance as a 
function of the Fermi wave vector for the unperturbed struc- 
ture (sohd line) and the linear correction (dashed) for xq — 0, 
3/0 ~ 4. Bottom: Similar, but for d^g/dV^ on the first plateau 
(point p). Inset: Gray-scale plot of the same data as a func- 
tion of xo and yo and the configuration with a minimal width 
Ay = 20 used in the numerical calculations. Anderson units 
are used for xo, yo, Ay, fcp, and Vt. 



The second-order correction [Eq. (12)], dominant on the 
plateaus, is always negative, exhibits small fringes, and 
has an overall decay with Xq , in agreement with Ref . [5] . 

Despite the simplicity of the harmonic saddle potential 
and its absence of mode mixing, the results obtained are 
quite representative of a generic QPC's behavior. This 
universality stems from the above mentioned transfor- 
mation into the transmission eigenmode basis, and was 
checked through numerical simulations in a variety of 
confining potentials. The results from recursive Green 
function algorithms 151 for the double harmonic oscil- 
lator (not shown) reproduce our analytical calculations, 
while those corresponding to a tight-binding model of 
a hard wall QPC are presented in Fig. [l] The change 
dg/dV^ is shown in the upper panel as a function of the 
tip position for a Fermi energy in the first conductance 
step (point s), and as a function of kp in the inset. Con- 
firming our analytical prediction, it is significant only 
in the step regions. The lower panel shows the change 
d^g/dV^ on the first plateau (point p). In agreement 
with Eq. (10), it is a negative correction that decreases 



with Xq and j/o further from the constriction center. The 



inset shows the full dependence of d^g/dV^ on the tip 
position. A checkerboard pattern consistent with exper- 
imental findings [4] can be observed in the second-order 
correction on the first plateau. The choice of a ^ tip 
was made for simplicity, but a smoother tip can be easily 
implemented. 

SGM measurements can be considerably influenced 
by disorder near the QPC. The conductance changes 
[Eqs. ([9]) and ( 10 )] fully incorporate the disorder through 
its effect on the scattering amplitudes and the scattering 
wave-function basis for expressing V. The transformation 
in mode space that singles out the propagating modes 
illustrates that the conclusions regarding the relative im- 
portance of g^^' and g"' remain valid in the presence 
of disorder. Since a disordered wire containing a QPC 
can be mapped onto an unconstrained conductor with 
a renormalized mean free path |27| . the random-matrix 
theory of quantum transport [21] can be used to deter- 
mine the statistical properties of the SGM patterns far 
from the QPC. The use of semiclassical methods 20 to 
evaluate the conductance corrections [Eqs. (|9| and ( [To| ] 
can yield an intuitive interpretation in terms of classi- 
cal trajectories bridging the gap between classical [U [5] 
and quantum [THIO] calculations and providing an un- 
derstanding of the region near the QPC (where coherent 
branching and transient behaviors occur). 

In summary, expressions for the lowest-order correc- 
tions to a nanostructure's conductance caused by a weak 
perturbation have been given. Applied to the SGM study 
of a QPC, the first-order corrections in the tip strength 
dominate at the conductance steps and are suppressed on 
the plateaus. The second-order corrections become dom- 
inant on the plateaus. The two expressions are of great 
interest since, even if most of the existing data have been 
obtained in the plateau regions, recent experiments 0] 
also explored the behavior at the mode opening finding 
a different interference pattern. It is also in this regime 
where the interaction effects [IHlIin] are expected to show 
their signature. 

The currently used tip potentials are strong enough to 
create a divot of depletion in the two-dimensional elec- 
tron gas [IHlj, which is clearly beyond the perturbative 
regime. However, we anticipate that much of the in- 
terpretation and many of the lessons derived from the 
perturbation calculations would still apply. In our quest 
for a rigorous interpretation of the SGM measurements, 
the perturbative approach is an unavoidable landmark 
towards the understanding of this fascinating problem. 
Note that in experiments on nanostructures which do not 
present conductance quantization, and for which there- 
fore the linear correction [Eq. ^] dominates the SGM 
response, the tip voltage has been varied over large in- 
tervals [SI [Zj exhibiting the expected linear behavior. 
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